The Double Life of Thermal QCD 



Mohammed Mia, Keshav Dasgupta, Charles Gale, and Sangyong Jeon 

Department of Physics, McGill University, 3600 rue University, Montreal, QC, Canada H3A 2T8 



Abstract 

We study the gravity dual of a thermal gauge theory whose behavior parallels that of thermal 
QCD in the far IR. The UV of our theory has infinite degrees of freedom. We holographically 
renormalize the supergravity action to compute the stress tensor of the dual gauge theory with 
fundamental flavors incorporating the logarithmic running of the gauge coupling. From the 
stress tensor we obtain the shear viscosity and the entropy of the medium at a temperature 
T, and investigate the violation of the bound for the viscosity to the entropy ratio. This is a 
shortened and simplified companion paper to hep-th/0902.1540, and is based on the talks given 
by M. Mia at the "Strong and Electroweak Matter 2008" workshop and the McGill workshop 
on "AdS/CFT, Condensed Matter and QCD" in the fall of 2008. 
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1. Introduction 

Strongly coupled Quark Gluon Plasma (QGP) poses theoretically challenging yet ex- 
perimentally accessible questions. The formation of QGP at RHIC is an example where 
theoretical descriptions are completely lacking at low energies because our perturbative 
techniques fail at strong couplings. However probing the nonperturbative regime of large 
N gauge theory through a gravity dual has led to some interesting results for the physics 
of the quark gluon plasma. A popular approach so far has been the AdS/CFT corre- 
spondence pQ, even though QCD is not a conformal field theory at UV. However, for 
certain gauge theories with running couplings, there exist gravity duals. At zero tem- 
peratures these gravity duals are studied with [2] and without [314] fundamental flavors. 
On the other hand, at high temperatures, there are examples of gravity duals without 
fundamental flavors in the literature [7j. In this paper we give an example of a gravity 
dual for a specific thermal gauge theory with fundamental flavors and logarithmically 
running coupling constants (see also [8]). Our aim here is to study quantities such as 
shear viscosity, entropy and the viscosity to entropy ratio with our gravity dual. 
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2. The Duality 



The gauge theory that we consider arises from stacking N D3 branes at the tip of a six 
dimensional conifold with base T 1,1 = S 2 x S 3 , then placing M D5 branes in such a way 
that it wraps the two cycle S 2 on the base of the conifold alongwith Nf D7 branes em- 
bedded via so-called Ouyang embedding [2] to introduce fundamental matter [5] . Finally 
by euclideanising and periodically identifying the time coordinate we can introduce tem- 
perature in the gauge theory. However the high temperature, i.e temperature above the 
dcconfincment temperature that we would be interested in, will be generated by inserting 
a black hole in the dual gravitational backgound [5]. The D3 brane world volume and 
the unwrapped directions of D5 and D7 branes extends in the four Minkowski directions. 
Strings can end on any of these branes and the excitation of the strings are described by 
an SU(N+M) x SU(N) gauge group with Nf flavors. If gi, gi denote the gauge couplings 
of SU (N + M) and SU (N) , then they have nontrivial beta functions (at zero temper- 

at-e) ^ + £ ] = and ^ - = 3M (l + lo gA ) . We 

see that the two gauge couplings run in opposite directions and in the IR SU(N + M) 
flows to strong coupling. Performing a Seiberg duality transformation, we identify the 
strongly coupled SU(N + M) with a weakly coupled SU(N - (M - Nf)) at the IR. We 
see that not only the number of colors are reduced, but the difference of the size of the 
gauge group now decreases from M in [3] to M — Nf. This difference will decrease by the 
increments of Nf until it is smaller than or equal to Nf. Then there are two possible end 
points: (a) if N is still greater than zero then we will have an approximately conformal 
theory, or (b) if N decreases to zero but with finite M left over then we will have a 
SU(M) theory with Nf flavors that confines in the far IR (see [5] for more details). The 
latter theory, or more particularly the high temperature limit of the latter theory, is what 
we are interested in and henceforth we will only consider thaiPI. 

At zero temperature the dual description (assuming of course that in this limit the 
gauge theory decouples from gravity) involves D7 branes, deformed conifold and fluxes 
|2J. Due to the underlying RG flows the supergravity description does not capture the 
choppy cascading nature of the gauge theories! Rather the supergravity description only 
captures the smooth RG flows in the theory. This would also mean that the actual number 
of colors in the theory is rather subtle to define. These details are explored in |4|5j . 

Once we switch on a temperature in the gauge theory, the dual gravity description 
loses its simplicity. We can no longer claim that the fluxes, warp factor etc would remain 
unchanged. Even the internal manifold cannot remain a simple warped deformed conifold 
any more. All the internal spheres would get squashed, and at r = there could be both 
resolution as well as deformation of the two and three cycles respectively. Due to this 
complicated nature of our background our entire analysis in [5] is based on the following 
limit: 

(g„ g s N f , glMNf, 9 -^-\ -> 0, (g s N, g s M) -> oo (1) 



1 For a review of Seiberg dualities and cascading theories, see [4]. For a review on brane constructions 
for cascading theories see [B]. 
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In [5] we presented our results to 0(g s Nf, g s M 2 /N), and discuss how to extend this 
to higher orders in g s Nf and g s M 2 /N. In the limit where the deformation parameter is 
small, we show that to 0(g s Nf, g s M 2 /N) we can analytically derive the background tak- 
ing a resolved conifold geometry. The resolution parameter a depends on g s Nf,g s M 2 /N 
as well as on the horizon radius i.e a = do + 0(g s Nf, g s M 2 /N,rh) with ao constant. 
Our background then has the following form: 



ds = 



1 

Vh 



F 3 = 2M Ai (^1 

3g s MN f 

in 

3g s MN f 

8tt 



- gi{r)dt 2 + dx 2 + dy 2 + dz 2 



3g s N f 



Vh 



g 2 (r)- i dr 2 + dM 2 5 



2n 
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where T>(r,ip) = — A e^ — \de^\ gi are the black hole factora[f]; e a = 3s ^ f \ Aj,B.j are 
C( a2 ; 9sNf) corrections that take us away from the zero temperature Ouyang background 
[2], and are worked out in details in [5]; and is a warped resolved conifold with the 
warp factor: 



h = 



L 4 



2L 4 
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where ej,rv a ) etc are defined as: 
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and I/ 4 is defined in the same way as in the AdS/CFT correspondence. The D7 branes 
are embedded via Ouyang embedding in this background. Note that our analysis to this 



Note that they are generically unequal. 
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order is analytic, but beyond this order we loose all analytic control and we can only find 
the background numerically (see also [7j for the background without flavor D7 branes 
and [9] for the background with flavor branes) . Furthermore the log r dependences of all 
the fluxes, warp factor and dilaton is only for regions close to the D7 brane. To get a 
good behavior at far infinity we need to embed the whole system in F-theory [TDJ . That 
would mean inserting extra 24 — Nf seven branes at infinity; and then all the fluxes 
would go as inverse powers of r at large r. In fact in [5] we made a concrete prediction 
for the large r warp factor: it is given by ([3]) except that the sum over a now ranges from 
— oo to +oo. This way we can get rid of possible Landau poles (or naked singularities) in 
the background (see also [H]). On the other hand for finite g s Nf,g s M 2 /N the small r 
behavior is almost identical to (J3j) . However we couldn't derive the precise form for the 
large r behavior of the warp factor. Therefore in [5] we only speculated a possible inverse 
r behavior for the warp factor. 

Once we know the dual background we can say that the Hilbert space of the gauge 
theory can be obtained from the Hilbert space of the string theory on this geometry. In 
fact the dual supergravity background captures the strongly coupled gauge theory where 
we expect the RG flow to be smooth without any choppy Seiberg dualities, although 
we might be able to interpret the strongly coupled gauge theory also as some kind of 
approximate cascading theory 

As we now know very well, the UV completion of cascading type theories require 
infinite degrees of freedom. Once we have infinite degrees of freedom at the UV, we no 
longer expect a finite boundary action from supergravity analysis! What we need is to 
regularise and renormalise the supergravity boundary action so that finite correlation 
functions could be extracted. This would also mean that the usual Witten type proposal 
[12] for the AdS/CFT correspondence can be re-expressed in terms of the boundary 
variables to give us the complete picture. Therefore we can rewrite the ansatze proposed 
by Witten et al 12J for our background geometry to take the following Wilsonian form: 



where M is Minkowski manifold, Stotai is the low energy type IIB action defined in 
the string framePI. <f>o is the fluctuation over a given background, Squ is the Gibbons- 
Hawking boundary term [14] and S'counterterm is the counter-term action added to renor- 
malise the action. 

The interesting part now is that we can define two classes of theories in this kind of 
backgrounds by introducing a cut-off at r = r c : 

• The first class is to analyse the theory right at the usual boundary where r c — ► oo. 
This is the standard picture where there are infinite degrees of freedom at the boundary, 
and the theory has a smooth RG flow from UV to IR till it confines (at least from the 
weakly coupled gravity dual). 

• The second class is to analyse theories by specifying the degrees of freedom at generic 
r c and then defining the theories at the boundary. All these theories would meet the 
cascading theory at certain scales under RG flows. The gravity duals of these theories 

3 There are subtleties involved in using •S*° r t i ^ 1 1 g ~ ■S'guQj^A' ^ or detailed discussions, consult [5]. 





exp (Stotai [</>o] + SqH + •S'counterterm) 



exp S re n[<M 



(5) 
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are the usual resolved-deformed conifold geometries cut-off at various r c with appropriate 
UV caps added (or, alternatively, degrees of freedom specified). These UV caps are non- 
trivial geometries added to the resolved-deformed conifold geometry from r = r c to 
r = od 4 1. In addition to that, both these classes of theories ought to have appropriate 
number of seven branes so that all the Landau poles could be removed. 

We can also say that the boundary degrees of freedom are now given by Af uv where 
■A/uv = e~",e — > and n >> 1 specifies the parent cascading theory. For all other 
theories that we define at the boundary r = oo, by cutting off the geometry at r = r c , 
will have Af uv going as jV uv = e~" but now n > 1. Yet another interesting thing about 
all these theories is that any thermodynamical quantitities defined in these theories will 
be completely independent of the cut-off r c and only depend on the temperature T (with 
an exponentially small dependences on the boundary degrees of freedom)! This will be 
briefly demonstrated in the next section (see [5] for the full exposition). 



Fig. 1. A very approximate way to denote all the UV completed theories that we could study in our 
framework by cutting off the geometry at various r = r c . For a more correct illustration of the situation 
see [5]. The outermost blue curves on both sides denote the parent cascading theory. 



3. The Stress Tensor 



One of the important thermodynamical quantity to extract from our background would 
be the stress tensor. As the stress tensor couples to the metric, it's expectation value is 
obtained from the renormalised action S Ten via (T i] (A c )) = S -^M w here kj is the 

four dimensional effective metric induced by the ten dimensional metric of the form @ by 
approximating the metric first by a five dimensional slice and then taking the required 



4 For example we could add AdS geometries from r = r c to r = oo. This is like adding Af = 4 or 2 degrees 
of freedom at the cut-off scale to UV complete the gauge theories. One issue here is the connection to 
the work of |13| . The question is can we write a boundary theory at r = r c itself instead of going to the 
actual r = oo boundary? As shown in 1 131 this is possible in AdS case because the theory on the surface 
of a ball in AdS space doesn't have to be a local quantum field theory. The non-local behavior in such a 
"boundary" theory is completely captured by the Wilsonian effective action at the so called boundary. 
As detailed in [5], this is tricky in the Klebanov-Strassler model precisely because there is no unique 
strongly coupled gauge theory here. Again, as in the RG flow pictures of [5], when the dual gravity 
theory is weakly coupled we have a smooth RG flow in the gauge theory side, but the theory at any 
given scale can be given by infinite number of representative gauge theories none of which completely 
capture the full dynamics. Thus it makes more sense to define the theories at r = oo boundary only and 
not at any generic r = r c . 
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pull-back (see [5] for details on the procedure). At every fixed r the five dimensional 
metric induces a four dimensional metric and with appropriate rescaling of the time 
coordinate, we can view this effective four dimensional metric as a flat Minkowski metric 
(with possible perturbations). If 0^ denote the perturbations on the background by 
quark strings that stretch from the D7 branes to the horizon of the black hole, then 
writing the supergravity action up to quadratic order in Oij ; using integration by parts 
together with appropriate Gibbons-Hawking terms; and holographically renormalising 
the subsequent action, we obtain: 



T 



a (4)L8] 



_ f d q sr-y I tjrmn , fjnm) „(4)[/3] _ A(KV l ? 4- K7"?W 
im+quark — J (2^ d ' ' ^ ^ lctl \ a \ ' nn * l«l M ' T 

oo ^ 

+ (K^ + Kffi)sMW + b^ 3) J n 5 nm e-^ + 0(Te-^) (6) 

3=0 J 

where (b^yAf uv ) together will specify the full boundary theory for a specific UV com- 
plete theory. The quantities }j mn , J m ,K mn etc. with m,n = 0,..,3 are completely in- 
dependent of the radial coordinate r and depend on the 3+1 dimensional spacetime 
coordinates. Note also that the stress tensor depends only on the temperature T and is 
independent of any cut-offs. The dependence on the UV degrees of freedom is exponen- 
tially suppressed, and in the limit Af uv — e~ n ,n » 1 we reproduce the result for the 
parent cascading theory. 



4. Results and Discussion 



With the general formulation of stress tensor ([6]), which is similar to the AdS/CFT 
results 16 1 in the limiting case M — Nf — 0, we can compute the wake a moving quark 
creates in a plasma with Oij being the metric perturbation due to string (see [17] for an 
equivalent AdS calculation) . Furthermore we can compute the shear viscosity rj from the 
Kubo formula with the propagator obtained from the dual partition function ^ and the 
entropy density s using Wald's formula. The result for the ratio rj/s is given by: 

s 4 7 r+Mog 2 (l-T 4 e- 4A/ '-) 

B (4tt 2 - log 2 C ) + 4ttA log C D 
(47r 2 -log 2 C ) 2 + 16^ 2 log 2 C 

where (A , B a , C D , a^j are constants that depend on the temperature T and e _JV " v ; and 
c 3 is the coefficient of the Riemann square term coming from the backreactions of the 
D7 branes in the background. These have been explicitly worked out in 15]. The C3 
dependence of the rj/s ratio first appeared in [18] . 

With certain choices of parameters (g s = 0.01, N = 1000, M = 100, Nf = 3, c 3 = 
0.001) and ignoring 0(g s Nf, g s M 2 /N) corrections to the black hole parameters <?j(r) the 
result for rj/s is shown in figure 2. It is easy to see that we are violating the celebrated 



C3K 

3L 2 (1 - T 4 e~ 4A H 3 
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KSS bound [TS0. The plot in figure 2 is shown for r\j s with and without Riemann Square 
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Fig. 2. Plot of rj/s versus the UV degrees of freedom. 

term. The x-axis is defined as k — where AT UV is the UV degrees of freedom and T 
is the temperature of the cascading theories. For the parent cascading theory, and their 
corresponding family of theories, k — > oo and we see a violation of the bound (the solid 
blue line). As k decreases (assuming this is possible!) the red dashed line dips slightly 
below the 1/Att axis, but the solid blue line remains considerably below the 1/4-/T axis. 
For k sufficiently small the bound is not violated. However all the models that we studied 
here and in [5] can only realise the infinite k limit, so that the small k limit depicted in 
Fig. 2 above is just an extrapolation of ([7]). 



5. Conclusion 

We have analyzed the gravity dual of a thermal gauge theory with logarithmic running 
coupling which may resemble QCD in the far IR. It appears that the KSS bound for 
r)/s may be violated in certain gauge theories if one carefully takes into account both 
the running of the couplings as well as the backreactions of the flavor branes in the 
dual gravity picture, for a range of parameter values (see details in [5]). It is of course 
important to test and verify the robustness of this limit, and we view the current work as 
contributing to this effort. More work is needed in order to identify the size and extent of 
this violation, at the moment this violation is only parametric and our parameters need 
to have better defined physical origins. In the end, one may need to rely on the empirical 
identification of key quantities, like transport coefficients for example |21j . In this regard, 
the role of heavy ion experiments at RHIC and at the LHC can't be overestimated. 



5 Recent papers dealing with the violation of the KSS bound are 18 20 . 
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